Noise Intensity-Intensity Correlations and the Fourth Cumulant of Current 

Fluctuations 
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We report measurements of the correlation between intensities of noise at different frequencies 
on a tunnel junction under ac excitation, which corresponds to two-mode intensity correlations in 
optics. We observe positive correlation, i.e. photon bunching, which exist only for certain relations 
between the excitation frequency and the two detection frequencies, depending on the dc bias of the 
sample. We demonstrate that these correlations are proportional to the fourth cumulant of current 
fluctuations. 



Intensity-intensity correlation is a very powerful expe- 
rimental tool which allows the characterization of statis- 
tical properties of light ,1]. In particular, it is the quan- 
tity that differentiates classical- from quantum light, by 
showing evidence of photon bunching or anti-bunching. 
Historically, the first experimental setup which addressed 
such correlations is that of Hanbury-Brown and Twiss Q 
on the broad spectrum of a distant star, whose light was 
measured with two spatially separated detectors. This 
setup has been repeated many times to characterize light 
sources (often with narrow spectra) using a beam splitter 
and again two detectors. Powerful theoretical tools have 
been developed to understand these measurements, in 
particular the famous second order correlation function 
(/ {t) I {t + t)), which precisely characterizes the correla- 
tion between the intensity I{t) of light at time t with that 
of light at time t + t. This quantity can be expressed in 
terms of the creation/annihilation operators of the pho- 
ton field, taking into account the type of detector that is 
used (photodetector, homodyne detector, etc.) Q. 

In mesoscopic physics, a lot of effort has been put in the 
measurement and understanding of current fluctuations 
0, [3| • Current or voltage fluctuations are simply another 
point of view of a fiuctuating electromagnetic field, which 
in the optics community are rather thought of as time 
dependent electric and magnetic fields. Thus, a noisy 
electronic device is also a light source, and it is natural to 
apply tools developed in optics to analyze high frequency 
electronic noise [6|-|8|. 

In the following we present a measurement of the cor- 
relation G2 = {P1P2) — (Pi) {P2) between the powers Pi 
and P2 of light at two different frequencies in the GHz 
range, /i and /2, as proposed in [8]. The source of mi- 
crowave light is a tunnel junction that is dc biased and 
excited at frequency fo. We observe that Pi and P2 are 
correlated when /q = /i + /2, /2 - /i or (/2 - /i)/2, de- 
pending on the dc bias voltage of the sample. We demon- 
strate that G2 is simply proportional to the fourth cumu- 
lant of the current fluctuations generated by the tunnel 
junction. We have calculated this quantity for various 
excitation voltages and found very good agreement with 
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FIG. 1. Experimental setup. Diode symbols represent fast 
power detectors. 



our data. 

Experimental setup. We have chosen to perform the 
measurement on the simplest system that exhibits well 
understood shot noise, the tunnel junction. The sample 
is a i? = 22 il Al/Al oxide/ Al tunnel junction, similar to 
that used for noise thermometry [ol, cooled at T = 3.0 K 
so the aluminum remains a normal metal. The experi- 
mental setup is depicted on Fig. [TJ The combination of 
a triplexer and a bias tee allows the junction to be con- 
nected to different instruments depending on frequency. 
The sample is voltage biased by a dc source and con- 
nected to two microwave signal generators, one ranging 
from 10 MHz to 4 GHz and the other above 8 GHz. In 
the 4 — 8 GHz frequency range, the noise generated by 
the junction at point A is amplified by a cryogenic am- 
plifier. The setup thus allows the noise to be measured in 
the 4 — 8 GHz range while ac excitation can be achieved 
both below 4 GHz and above 8 GHz. This insures that 
the amplifier never sees the ac excitation of the sample, 
which avoids spurious signals due to non-linearities in the 
amplifier at high excitation power. 
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FIG. 2. (color online) Reduced (32 vs excitation frequency /o 
for two bias voltages: Vdc = (in red) and Vdc = 2.24 mV (in 
blue). The ac excitation amplitude is Vkc = 1.1 mV for both 
measurements. The exact shape of the peaks is determined 
by the frequency response of the experimental setup. 



The amplified noise is split into two branches: in 
branch 1 (respectively branch 2) the signal is bandpass 
filtered around /i = 4.50 GHz with a bandwidth A/i = 

0. 72 GHz (rcsp. /a = 7.15 GHz, A/2 = 0.60 GHz). In 
each branch a fast power detector (diode symbol on Fig. 
[1]) measures the "instantaneous" (averaged over a few ns) 
noise power integrated in the corresponding bandwidth, 

1. e. VBf.{t) oc Pk{t) with fc = 1,2. The dc voltage at 
point Bk is thus proportional to the noise spectral den- 
sity, S{fk) = (|*(/fc)P) with i{f) the Fourier component 
of the instantaneous current: {vb^) oc S{fk)^fk- A dc 
block (capacitor on Fig. [T|) removes the dc part of to 
give vc^{t) oc Pk{t) - (Pfc) = SPk{t). Finally, voltages 
at points Ci and C2 are simultaneously digitized with 
a 2-channel, 14-bit, 400 MS/s acquisition card and the 
correlator G2 = {SP1SP2) is computed in real time by 
a 12-core parallel computer, as are the autocorrelations 
of both channels ^JP^). The gain has been adjusted to 
draw on the full dynamic range of the digitizers. 

In order to have a quantitative measurement of G2, it 
is necessary to calibrate the ac excitation voltage reach- 
ing the sample and the overall gain of the detection. The 
calibration of the ac voltage across the sample is per- 
formed by measuring the usual photo-assisted noise, i.e. 
S vs Vdc, in the presence of a microwave excitation for 
various excitation voltages [ll| . The temperature is 
large enough (fc^T ^ hf 1^2) to approximate the noise 
measured in either branch by its value at zero frequency. 
In the absence of ac excitation, the noise spectral den- 
sity is given by S'o(Vdc) = eGVdc coth(eVdc/2fcsr) with 
G = 1/R, the sample's conductance. Since the excitation 
frequency /o is always such that hfo <C fcsT, the photo- 
assisted noise can be approximated by its time-average 
value as if the junction were responding instantaneously 
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FIG. 3. (color online) Reduced G2 vs dc bias voltage for var- 
ious ac excitation amplitudes at frequency fo ~ f+ = 11.6 
GHz. Symbols are experimental data and solid lines theoret- 
ical expectations of Eq. (HI). 

to the time-dependent voltage: 

^oC^dc + KcCOS0) — . (1) 

All measurements of S vs Vdc at fixed V"ac are indeed 
identical and behave according to Eq. ([T]) whether fo = 
1 kHz or /o = 12 GHz (data not shown). The setup gain 
calibration will be discussed later. 

Results. The first step for the study of G2(Vdc, V'ac, fo) 
is to determine the excitation frequency for which G2 ^ 
0. This is done by measuring G2 at fixed bias volt- 
age (Vdc = or Vdc = 2.4 mV) and fixed excita- 
tion amplitude (Vac ~ 1-1 mV) while varying the ex- 
citation frequency /q from 10 MHz to 3.83 GHz and 
from 10 GHz to 14 GHz. These results are reported in 
Fig. [21 where G2 has been reduced to units of as 
noise is often given in terms of equivalent temperature 
Tnoisc = RS/2kB. We observe that G2 is large at low 
frequency for both bias voltages, whereas at higher fre- 
quencies G2 strongly depends on Vdc. For Vdc = 2.4 mV, 
G2 shows peaks at /o ~ /_ = /2 — /i = 2.6 GHz and 
fo — f+ ~ f2 + fi ~ 11-6 GHz. At zero bias, G2 peaks 
at fo ~ /-/2 = 1.3 GHz. Choosing fo = /+/2 was not 
possible with this experimental setup since it lies in the 
4 — 8 GHz range. 

We then consider the variation of G2 as a function 
of the dc voltage for various ac excitation amplitudes 
at fixed excitation frequencies, chosen according to the 
position of the peaks we observed. Data on Fig. [3] corre- 
spond to an excitation at /o = /+. We observe that G2 
is maximal at high dc bias and vanishes at Vdc = 0. We 
obtained identical results for fo = /- (data not shown). 
Data on Fig. [4] correspond to fo — /-/2. Here G2 peaks 
at dc bias and decays when |Vdc| increases. Moreover, 
the maximum of G2(/o = /±) for a given Vac is an or- 
der of magnitude larger than that of G2(/o = /-/2). G2 
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FIG. 4. (color online) Reduced G2 vs dc bias voltage for var- 
ious ac excitation amplitudes at frequency /o = /-/2 = 1.3 
GHz. Symbols are experimental data and solid lines theoret- 
ical expectations of Eq. ((H). 



measured for low frequency excitation (data not shown) 
is similar to the behaviour at /o = /±, Fig- El for small 
Vac but differs slightly at higher excitation voltages, as 
will be discussed below. 

The fourth cumulant of noise. In order to describe 
and understand the meaning of the measured correla- 
tor, let us express the noise generated by the sample 
in Fourier space. Introducing the Fourier component 
of the power fluctuations SPk{e) at frequency e, one 
has: G2 — J de {SPi{e)SP2{—e)). Here the integral 
spans over the output bandwidth of the power detectors, 
Ae ~ 100 MHz, except for very low frequencies which 
are removed by the dc block. Specifically, one always 
has e ^ 0. Power fluctuations are related to current fluc- 
tuations by SPkie) <xjdf [iif)i{e - f) + i{-f)i{e + /)]. 
This integral spans over the bandwidth of the bandpass 
filters, i.e. fk ± A/fe/2. Thus, G2 is proportional to the 
correlator between currents at four different frequencies 
which, by definition, is the time averaged fourth cumu- 
lant of current fluctuations taken at frequencies /i, /2 
and e (with e ~ 0). 

G2 (X {t {fi)i (e - A) z (/2) t {~e ~ f2)) A/iA/2A£ (2) 

This represents the fourth cumulant {{i^{t))) = (*^(^)) — 
3 («^)^ of the current when it is filtered to keep only fre- 
quencies close to ±/i or ib/2 and averaged over time t. 
In order to detect the fourth cumulant and not the fourth 
moment of current fluctuations, it is crucial that the four 
frequencies be different so that correlators («(/)«(—/)) 
are never involved. Experimentally, the separation of 
the signal into two branches with non-overlapping band- 
widths insures that /i ^ if 2, while the presence of the 
dc blocks, by imposing e ^ 0, prevents all other possible 
occurrences of such correlators. 

Single channel autocorrelations (SPi) are 



also related to fourth order current correlators 
{i if) i{e-f)i if) I (-£ - /')). However, here / and /' 
belong to the same frequency band, so the correlator is 
dominated by terms / = — /', and (SP^) oc {Pkf' . The 
fourth cumulant G2 is only a very small correction to 
this. Thus, power correlations within the same frequency 
band are totally dominated by Gaussian fluctuations. 
Experimentally, we indeed observe that 
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(SP^) oc (Pk)'^, which provides a first way to calibrate 
our data. Note that since the amplifier noise dominates 
(Pk), it also dominates (SP^) . In contrast, only the 
fourth cumulant of the amplifier's noise contributes to 
G2. 

Interpretation. We will first focus on the low frequency 
excitation. When the voltage oscillates very slowly 
(/o < Ae), the outputs of the power detectors are pro- 
portional to the "instantaneous" noise spectral density 
S{t) = So (Vdc + Vac cos {27: fat)). If /o is large enough to 
pass through the dc blocks (/o > 100 kHz), G2 is given by 
the time average of 5S'^{t) with 5S{t) = S{t) - S. Here S 
is the time-averaged value of S{t), i.e. the photo-assisted 
noise of Eq. ([ij . In terms of the Fourier components of 
the current, the amplitude of the oscillation of 5S{t) at 
frequency n/o is given by X„(/, fo) = {i (/) i (nfo - /)), 
with n positive or negative integer. This correlator has 
been calculated and measured in the quantum regime, 
at very low temperature 13-1^. In the high tempera- 
ture, classical regime that corresponds to the present ex- 
periment, Xn reduces to the Fourier coefficient of SS{t) 
oscillating as exp{in0) with 9 = 27r/oi, i.e.: 



Xn — 



dO 

So {Vdc + Kc cos 9) exp (inO) — (3) 

ZTT 



which is independent of / and /q. The measured fourth 
cumulant is then given by G2 {fo ^ ^s) = K X^ 
which, according to the Parseval-Plancherel relation, is 
the time average of 6S{t)'^. The normalization factor 
K ~ 4A/1A/2 corresponds to the integral over the fre- 
quencies that contribute in the setup. For small Vac, 
6S{t) ~ Vac^^ cos(27r/ot) oscillates at frequency /o, but 
for large Vic, SS{t) may contain many harmonics of fo. 
Those harmonics are filtered out as fo approaches Ae, 
so G2 should vanish for fo » Ae, as observed in Fig. 
[21 These formulae reproduce our results very well (data 
not shown) and provide another way to calibrate the 
overall gain of our setup. This calibration and the one 
based on (<5P^) are consistent within 10%, even though 
they involve integrals over frequencies which slightly dif- 
fer. In particular, the phase difference between the two 
branches of the setup, which has been adjusted to max- 
imize the signal, does not infiuence the autocorrelations 
(SP^) while it affects the cross-correlation {5P1SP2). 

In terms of probability distribution of the current fluc- 
tuations one can, at low excitation frequency, think of 
the amplitude of the current fluctuations as being slowly 
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modulated by the ac voltage: if one considers a Gaus- 
sian with a variance oscillating in time, the time-averaged 
probability distribution is not a Gaussian and has a non- 
zero fourth cumulant. 

To understand why G2 may be nonzero at high ex- 
citation frequency, let us first discuss the correlator 
(*(/)*(/'))■ This will be nonzero only for frequencies 
such that f + f = nfo with n, an integer. The case 
n = corresponds to S whereas n describes the 
noise dynamics 



15j |. If we consider our experimental 
setup, relevant frequencies are close to ±/i and ±/2, so 
this condition becomes /i ± /2 = ri/o, i.e. fo — f±/n. 

In such cases, the fourth order correla- 
tor of Eq. ([2]) is dominated by the terms 

{^ih)^i±h)){^i-h)^{Tf2)) = Xl 
has: 



Thus, one 



G2{fo^ f±/n)^KXl 



(4) 



For small Kc, « i (^)' V^^ and X| « (^) 
This is why for fo = f± we observe on Fig. |3] that 
G2 — KXf is at Vdc = and maximal at large V, with 
G2(/o = /±,14c > fcsT/e) = X(eGKc/2)2, whereas for 
/o = /-/2, G2 — KX2 is maximal at Vdc = 0, with 
G2(/o = /-/2, T4c = 0) = A' {e'GVll2AkBT)\ and de- 
cays at finite dc bias, as seen on Fig. S) We also observed 
that G2(/o = /±) oc Kc and G2(/o = /-/2) oc Kt for 
small Kic (data not shown). Sofid lines on Fig. [3]and|4] 
represent the theoretical predictions of Eqs. ([3]) and (jl]). 

In order to make a connection with experiments per- 
formed in optics, let us define the dimensionless correla- 
tor: 
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(5) 



which is usually referred to as the same-time two-mode 
second order correlator of the electromagnetic field ra- 
diated by the junction. To calculate it, one has to sub- 
tract the contribution of the amplifiers from P^'. (Pk) oc 
(S+Samp) with Samp the Current noise spectral density of 
the amplifier. The result is plotted on Fig. [S]as a function 
of Vdc for Vac = 1.0 mV and both fo = /+ (red triangles) 
and fo = /-/2 (green circles). It is clear from these data 
that we always observe a positive correlation between 
the power fluctuations, i.e. photon bunching (172 > !)■ 
Each acquisition performed by the digitizer (integrated 
over T — 2.5 ns) corresponds to an averaged number of 
(n2) = P2T/hf2 ^ 21 photons at 7.2 GHz emitted by the 
sample at Vdc — Vac — 1 niV and fo — /+, plus ~ 40 
photons from the amplifier. Thus our experiment does 
not measure correlations at the single photon level, but is 
not very far from that limit, which can be reached by low- 
ering the temperature and the ac power. Note that our 
power detector, which measures the square of the elec- 
tric field, cannot differentiate absorption from emission 
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FIG. 5. (color online) Measured g2 vs dc bias voltage at V^c = 
1.0 mV for an excitation frequency fo = /+ (red triangles) 
or /o = /-/2 (green circles). The dashed line at 32 = 1 
represents independent photons, i.e. chaotic light, while 52 > 
1 corresponds to photon bunching. 



of photons by the sample. This is extremely important 
when comparing to theories such as 0^ which consider 
photon detectors. In particular, the correlations between 
jhoton detectors will involve other current correlators 
3]. Whereas G2 at excitation frequency /o = /i + /2 
probably corresponds to absorption of one photon of fre- 
quency fo and emission of two correlated photons, one 
at /i and one at /2, the case of excitation at frequency 
/o = /2 — fi probably corresponds to two photons being 
absorbed, one at fo and one at /i while one photon at 
frequency /2 is emitted. We indeed observe no difference 
in G2 for fo = f+ and /o = /-• 

Perspectives. Current fiuctuations generated by a tun- 
nel junction are known to be non-Gaussian, and thus 
have a non-zero fourth cumulant even in the absence of 
ac excitation [13] ■ This should contribute to our measure- 
ment as ((*^)) c>c e^S'o. In the absence of excitation, we 
indeed observe a non-zero, voltage-dependent G2. How- 
ever, the signal we observe is much larger than ((j^)), and 
we believe it to be an artefact of our experimental setup 
(slightly overlapping bandpass filters, cross-talk in the 
digitizer). This background signal is negligible at large 
Vac but has to be subtracted from the measurement at 
low Vac for the results to be well fitted by Eqs. ([3]) and 
(|4|). Measurements of non-Gaussian noise are also known 
to be affected by environmental contributions l^- 2^ • In 
[lit I these contributions have been explicitly exhibited by 
exciting the sample with microwave radiation. Similarly, 
our results allow us to estimate the environmental con- 
tributions to an eventual measurement of the intrinsic 
fourth cumulant. Each Fourier component of the fluc- 
tuating environment will contribute as ~ X^. At large 
vohage, this gives {SV^)e'^G'^AfiAf2 where {SV"^) is 
the total (integrated over a potentially huge bandwidth, 
much wider than the detection one) voltage noise experi- 
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enced by the sample. It follows that the intrinsic fourth 
cumulant may be dominated by the environmental con- 
tributions. 

Another interesting issue raised by our experiment con- 
cerns the physical interpretation of G2 in terms of tem- 
perature fluctuations. Since the noise can be expressed in 
terms of a temperature, noise correlations are equivalent 
to temperature fluctuations which, once integrated over 
the detection bandwidth of the correlator, have units of 
K^, as in Figs. 2-4. Let us suppose we measure G2 on 
a macroscopic resistor. Then both Pi and P2 would be 
proportional to the the same sample's temperature, so 
G2 would be a measure of the variance of temperature 
fluctuations, G2 oc {ST"^). These are non-zero for any 
sample of finite size. Thus the thermal noise of a resis- 
tor is gaussian only if the sample is infinite, such that its 
temperature does not fluctuate at all. In other terms, the 
blackbody radiation should exhibit correlations similar 
to the one we have observed, which implies that thermal 
light is not exactly chaotic. 

Conclusion. We have shown that the noise power of a 
photo-excited tunnel junction measured at two different 
frequencies is substantially correlated. This is crucial in 
demonstrating the non-gaussian nature of noise signals 
and a first step towards the study of the intrinsic fourth 
cumulant of current fluctuations. In terms of photons, we 
have shown that photo-assisted shot noise corresponds to 
photon bunching. 
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